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ABSTRACT 



We present the estimates of the order 0(af) QCD corrections to R(s) = cr toi (e + e~ — > hadrons)/ 
cr(e + e~ — > R T = T(r — > z/ T + hadrons) /Y{r — > u T V e e) and to the deep inelastic scattering 

sum rules, namely to the non-polarized and polarized Bjorken sum rules and to the Gross- 
Llewellyn Smith sum rule. The estimates are obtained in the MS'-scheme using the principle 
of minimal sensitivity and the effective charges approach. 
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1 Introduction 



During the last few years, essential progress has been achieved in the area of the calculation of 
the next-next-to-leading order (NNLO) QCD corrections to the number of physical quantities. 
Indeed, the complete NNLO 0(a^) QCD corrections are known at present for the character- 
istics of e + e~ — > hadrons process U, 0, r — »■ v r + hadrons decay ||, [|J and for the deep 
inelastic scattering sum rules, namely the non-polarized Bjorken sum rule (BjnSR) 0, the 
Gross-Llewellyn Smith sum rule (GLSSR) and the polarized Bjorken sum rule (BjpSR) 0. 
Amongst the physical information provided by these results is the estimate of the theoretical 
uncertainties of the corresponding next-to-leading-order (NLO) perturbative QCD predictions 
for R(s) = o~ to t(e + e~ — > hadrons)/a(e + e~ — > 0, T(Z° — ► hadrons) (see. e.g., 0), 

R T = r(r -f u T + hadrons)/r(r -> v T V e e) @, BjnSR §, GLSSR and BjpSR ||. 



In view of the fact that the precision of the experimental data for a to t{e + e~ — > hadrons), 
T(Z° — > hadrons) and for the structure functions of the deep inelastic scattering is continuously 
increasing, the problem of estimating the effects of the next-after-next-next-to-leading order 
(NANNLO) O(of) corrections to the measurable physical quantities now arises. 

The first attempt to build the bridge between the results of the concrete calculations of 
the order O(ct^) corrections to R(s) and the higher-order coefficients of the corresponding 
perturbative QCD series was made in Ref. JTT]. However, it was later demonstrated that the 
considerations of Ref. JTTJ have definite drawbacks fL2| , ||13|| . A certain step in the direction of 



more substantiated estimates of the order 0{pQ) QCD corrections was made in the case of R T 



in Ref. 0. These estimates are based on the tendency, observed in Ref. |L4[], of the scheme- 
dependent uncertainties of the perturbative QCD predictions for R(s) and R T to decrease as a 
result of taking into account the order 0(al)-terms. The foundations of Ref. [14] were further 



generalized in the process of phenomenological studies of the QCD predictions for R T , taking 
into account the concrete pieces of the higher-order perturbative QCD corrections [15j (for the 
earlier development of a similar technique, see Ref. 0]). 

It should be stressed that the analysis of Ref. does not allow us to fix the sign of the 
order O(a^) contribution to R T . In this work we generalize the QED ideas of Refs. [T7|], |T8| ] 
and [19[ to the QCD case and estimate the order 0(a^) and 0(af) corrections to R(s), R T 



and deep inelastic scattering sum rules with the help of the principle of minimal sensitivity 
(PMS) ]I7J and the effective charges (ECH) approach [p0| , which is equivalent a posteriori to 
the scheme-invariant perturbation theory Contrary to the conclusions of Ref. we 

argue that the application of these approaches allows one to control the value and fix the sign 
of the higher-order QCD corrections to physical quantities, once the preceding ones are known 
in the particular scheme. 
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2 The Description of the Formalism 



Consider first the order 0(a ) approximation of a physical quantity 



N-l 



D N = d a(l + J2 ^ (!) 



i-1 



with a = a s /7r being the solution of the corresponding renormalization group equation for the 
/3-function which is defined as 



da N ~ x 
V 2 j^ = (3(a) = -f3oa 2 (l+Y: *<*) • (2) 



In the process of the concrete calculations of the coefficients di,i > 1 and Cj,i > 2, the MS 1 
scheme is commonly used. However, this scheme is not the unique prescription for fixing the RS 
ambiguities. In both phenomenological and theoretical studies other methods are also widely 
applied. 

The PMS [[HJ and ECH [^] prescriptions stand out from these other methods. Indeed, 
they are based on the conceptions of the scheme-invariant quantities, which are defined as the 
combinations of the scheme- dependent coefficients in Eqs. (1) and (2). Both these methods can 
claim to be the "optimal" prescriptions, in the sense that they provide better convergence of 
the corresponding approximations in the non-asymptotic regime, and thus allow an estimation 
of the uncertainties of the perturbative series in the definite order of perturbation theory. 
Therefore, the applications of the "optimal" methods allow one to estimate the effects of the 
order 0(a iV+1 )-corrections starting from the approximations D c ^ t {a opt ) calculated in a certain 
"optimal" approach |[7j - |19j . This idea is closely related to the QED technique of Ref. |23j , 



which was used to predict the renormalization-group controllable /n(m /1 /m e )-terms in the series 
for (g — 2) M from the expression of (g — 2) M through the effective coupling constant d(m^/m e ). 
In our work we are using a similar technique to estimate the constant terms of the higher-order 
corrections in QCD. 

Let us re-expand D < ^ t (a op t) in terms of the coupling constant a of the particular scheme 

Df{a opt ) = D N (a) + 5Dfa N+1 (3) 

where 

5Df = n N (d t , a) - n N (df, cf) (4) 

are the numbers which simulate the coefficients of the order 0(a N+1 )-corrections to the physical 
quantity, calculated in the particular initial scheme, say the MS'-scheme. The coefficients VIn 
can be obtained from the following system of equations: 

^(D N + n N a N+1 )=0(a N + 2 ), 
^-(D N + n N a N+1 ) = 0(a N+2 ), i>2 (5) 



2 The application of this similar technique to the analysis of the five-loop approximations of (g — 2)^ and 
(g — 2) e will be considered elsewhere. 
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where the parameter r = /3o£n(fi 2 / A 2 ) represents the freedom in the choice of the renormal- 
ization point fi. The conventional scale parameter A will not explicitly appear in all our final 
formulas. 

The explicit form of the coefficients Q 2 and Q3 in which we will be interested can be obtained 



by the solution of the system of equations (5), following the lines of Ref. | 17fl . We present here 
only the final expressions: 

^2 = rf rfi(ci + di), (6) 

Q 3 = d di(c 2 + ici<ii - 2d\ + 3d 2 ) . (7) 

It should be stressed that in the ECH approach df CH = for all % > 2. Therefore one gets 
the following expressions for the NNLO and NANNLO corrections in Eq. (3): 

6Di CB = ih(d 1 ,c 1 ) (8) 

6D* CH = n 3 (d u d 2} c u c 2 ) (9) 
where fl 2 and fl 3 are defined in Eqs. (6) and (7) respectively. 

In order to find similar corrections to Eq. (3) in the N-th order of perturbation theory 
starting from the PMS approach |T7|], it is necessary to use the relations obtained in Ref. 



p4]l between the coefficients rf MS and cf MS (i > 1) in the expression for the order 0(dp MS ) 
approximation D^ MS (apMs) of the physical quantity under consideration. The corresponding 
corrections have the following form: 



5D^ MS = 5D* CH - ^ (10) 

5D™ S = SDi CH . (11) 
Notice the identical coincidence of the NANNLO corrections obtained starting from both the 



PMS and ECH approaches. A similar observation was made in Ref. JT8[ using different (but 
related) considerations. 



3 Input of the Analysis 

Consider first the familiar characteristic of the e + e~ — > 7 — ► hadrons process, namely the 
D-function defined in the Euclidean region: 

D <« = « 2 f (jf^ 8 (12) 

Its perturbative expansion has the following form: 

D(Q 2 ) = 3£Qj[l + a + d 1 a 2 + d 2 a 3 + d 3 a 4 + ...] 

+ (EQ / ) 2 [J 2 a 3 + 0(a 4 )] (13) 
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where Qf are the quark charges, and the structure proportional to (£Q/) 2 comes from the 
light-by-light-type diagrams. The coefficients d\ and d 2 , d 2 were calculated in the MS'-scheme 
in Refs. [7| and |l| respectively. They have the following numerical form: 

dfs « 1.986 - 0.115/ 

df 3 » 18.244 - 4.216/ + 0.086/ 2 , d 2 « -1.240 . (14) 



Following the proposals of Ref. [25|, we will treat the light-by-light-type term in Eq. (13) 
separately from the "main" structure of the .D-function, which is proportional to the quark- 
parton expression D QP {Q 2 ) = 3T,Q 2 j. In fact, one can hardly expect that it is possible to predict 
higher-order coefficients di, i > 3 of the second structure in Eq. (13) using the only explicitly- 
known term d 2 . Therefore we will neglect the light-by-light-type structure as a whole in all our 
further considerations. This approximation is supported by the relatively tiny contribution of 
the second structure of Eq. (13) to the final NNLO correction to the .D-function. 

The next important ingredient of our analysis is the QCD /3-function (2), which is known 
in the MS-like schemes at the NNLO level [p6~| |. Its corresponding coefficients read 



A> = (11 --/)-« 2.75 -0.167/ 

153-19/ 
Cl ~ 66-4/ 
m 77139- 15099/ + 325/ 2 
° 2 ~ 9504- 576/ • 

Using now the perturbative expression for the Z)-function, one can obtain the perturbtive 
expression for R(s), namely 

R(s) = 3ZQ 2 f {l + a s + r 1 a 2 s + r 2 a 3 s +r 3 a 4 s + ...} 

+ (EQ f ) 2 [f 2 a 3 s + ...] (16) 

where a s = a s /7r, and 

n = d x 

r - d ^ r - d 
r 2 — a 2 — , r 2 — a 2 

r 3 = d 3 - n 2 p 2 (dx + ~cx) . (17) 

The corresponding 7r 2 -terms come from the analytic continuation of the Euclidean result for 
the D-function to the physical region. The effects of the higher-order 7r 2 -terms were discussed 
in detail in Ref. |27[ . 



The perturbative expression for R T is defined as 

I'M? ds 



Rr = 2 1 W2 (l-s/M 2 ) 2 (l + 2s/M 2 )R(s) 
~ 3[1 + a T + rla 2 r + r T 2 al + rla A r + ..] 
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where a T = a s (M 2 )/ir and R(s) is R(s) with with / = 3, (HQ/) 2 = 0,3£Q 2 substituted for 
3S | V /f > | 2 and | Kd | 2 + | | 2 ~ 1. 



It can be shown that in the MS'-scheme the coefficients of the series (18) are related to 
those of the series (13) for the D-function by the following numerical relations (see Refs. Pf, 



1]): 



T \ = df s {f = 3) + 3.563 = 5.202 
r\ = df^if = 3) + 19.99 = 26.366 

rl = dJ3(f = 3) + 78.00 . (19) 

In order to estimate the values of the order 0(a 3 ) and 0(a 4 ) corrections to R(s) and R T , we will 
apply Eqs. (6) - (11) in the Euclidean region to the perturbative series for the D-function and 
then obtain the corresponding estimates using Eqs. (17) and (19). This is one of the origins of 
the disagreement of our considerations with those presented in Ref. 



We now recall the perturbative expression for the non-polarized Bjorken deep-inelastic scat- 
tering sum rule 



BjnSR = C Fl p - vp (x,Q 2 )dx 
J o 



2 

1 a(l + d 1 a + d 2 a 2 + d 3 a 3 + ...) (20) 

3 



where the coefficients d± and d 2 are known in the MS scheme from the results of calculations 
of Ref. and Ref. || respectively: 



df b w 5.75 - 0.444/ (21) 

dfs » 54.232 - 9.497/ + 0.239/ 2 . (22) 
The expression for the polarized Bjorken sum rule BjpSR has the following form: 



BjpSR = C gT- en (x,Q 2 )dx 
Jo 



-1 — 1 {l-a(l + dia + d 2 a 2 + d 3 a 3 + ...)} (23) 
3 g v 



where the coefficients di and d 2 were explicitly calculated in the MS scheme in Refs. [0] and 
respectively. The results of these calculations read 

df™ ~ 4.583 - 0.333/ (24) 

dfs ~ 41.440 - 7.607/ + 0.177/ 2 . (25) 

It should be stressed that since deep inelastic scattering sum rules are defined in the Euclidean 
region, we can directly apply to them the methods discussed in Section 2, and thus predict the 
values of the coefficients d% using Eqs. (9) and (7) without any additional modifications. 
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It is also worth emphasizing that, in spite of the identical coincidence of the NLO correction 
to the Gross-Llewellyn Smith sum rule 



GLSSR = - f 1 F£ p+Up (x, Q 2 )dx 
2 Jo 

= 3{1 - a(l + dia + d 2 a 2 + d 3 a 3 + ...)} (26) 

with the result of Eq. (24) [fLOfl , the corresponding NNLO correction differs from the result of 
Eq. (25) by the contributions of the light-by-light-type terms typical of the GLSSR 0: 

Wglssr = (d 2 )Bj P SR + d 2 , d~2 = -0.413/ (27) 

Since these light-by-light-type terms appear for the first time at the NNLO, it is impossible 
to predict the value of the light-by-light-type contribution at the NANNLO level using the 
corresponding NNLO terms as the input information. However, noticing that at the NNLO 
level the corresponding light-by-light-type contributions are small, we will assume that the 
similar contributions are small at the NANNLO level also. After this assumption our estimates 
of the NNLO and NANNLO corrections to the BjpSR can be considered also as the estimates 
of the corresponding corrections in the perturbative series for th GLSSR . Note, that the 
Pade predictions of the order 0(a A ) contributions to the GLSSR ||28||, which do not take into 



account the necessity of the careful considerations of the light-by-light-type terms, have definite 
drawbacks. However, certain other results of Refs. and |29| deserve comparison with the 
results of our studies. 



4 Outputs of the Analysis 

The estimates of the coefficients of the order 0(a 3 ) and 0(a 4 ) QCD corrections to the D- 
function, R(s), BjnSR and BjpSR/GLSSR obtained following the discussions of Section 2 with 
the help of the results summarized in Section 3, are presented in Tables 1-4 respectively. Due 
to the complicated /-dependence of the coefficients ^2,^3 in Eqs. (6) and (7), we are unable 
to predict the explicit /-dependence of the corresponding coefficients in the form respected 
by perturbation theory. The results are presented for the fixed number of quark flavours 
1 < / < 6 and are related to the MS scheme. The estimates of the NNLO corrections, 
obtained starting from both the ECH and PMS approaches, are compared with the results of 
the explicit calculations. 

Using the results of Table 1 and Eqs. (19), we are able to predict the NNLO coefficient of 

R T \ 

{rl)f CH a 25.6 (28) 
{rl)f MS » 24.8 . (29) 
One can see the agreement with the explicitly calculated result 

{rDws = 26 - 366 • (3°) 
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Considering this agreement as the additional a posteriori support of the methods used, we use 
the estimate of the NANNLO coefficient for the D-function with / = 3 numbers of flavours as 
presented in Table 1, namely: 

df = 27.46 (31) 
and predict the value of the NANNLO coefficient of R T : 

(r T 3 ) est « 105.5 (32) 



5 Discussions 



We are now ready to discuss the main results of our analysis. 

1. The predictions of the NNLO corrections to the D-function, BjnSR and BjpSR/GLSSR 
are in qualitative agreement with the results of the explicit calculations of Refs. |I| H and 
H and respect the tendency of the corresponding coefficients to decrease with increasing 
number of flavours. 

2. The best agreement of the NNLO estimates with the exact results is obtained for the case 
/ = 3. This fact, which we do not understand, supports the application of the method 
used for estimating the NNLO and NANNLO corrections to R T . 

3. The positive feature of the estimates as obtained by us is their scheme-dependence. Notice 
that since the methods used correctly reproduce the renormalization-group-controllable 



terms |21|, the transformation from the MS scheme to other variants of the MS-like 
scheme will not spoil the qualitative agreement with the results of the explicit calculations. 

The results of Table 2 demonstrate that the ir 2 effects give dominating contributions to 
the coefficients of R(s). 

From the phenomenological point of view, the most interesting results are the estimates 
of the NANNLO corrections to R T (see Eq. (32)) and R(s) for / = 5 numbers of flavours. 
Taking a s (Mz) ~ 0.12, we get the estimate of the corresponding NANNLO contribution 
to both T(Z° — > hadrons) and T(Z° — > qq): 

5Y z o « -97(a(M z )) 4 « -2 x 10" 4 . (33) 

It is of the order of magnitude of other corrections included in the current analysis of 
LEP data (see, e.g., 



6. Taking a s (M T ) « 0.36 [[|, we get the numerical estimate of the NANNLO contribution 
to R T : 

SR T w 105.5 a 4 « 1.8 x 10" 2 . (34) 



It is larger than the recently calculated power-suppressed perturbative |31|] and non 



perturbative 32 contributions to R T . 
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Our estimate for d^(f = 3) of Eq. (31) is more definite than the one presented in Ref. 
p9| , namely d 3 (f = 3) = 55^24, and than the bold guess estimate d 3 (f = 3) = ±25, given 
in Ref. Moreover, the related estimate of the NANNLO contribution to R T (see Eq. 



(32)) is more precise than those presented in Refs. || and fl5]], namely 5R T = ±130 a\ 



@ and 5R T = (78 ± 25)a 4 . fllS] , and is smaller than the result of applying the Pade 
resummation technique directly to R T , namely SR T = 133a 4 . [^9] . 



The qualitative agreement of the results of Tables 3 and 4 for BjnSR and BjpSR with the 



corresponding Pade estimates of Ref. [2~8~ can be considered as the argument in favour 
of the applicability of both theoretical methods in the Euclidean region for the concrete 
physical applications. Let us stress again that in the process of these applications the 
light-by-light-type structures, contributing say to R(s) and GLSSR, should be treated 
seperately. 



Notice also , that the application of the Pade resummation technique to R(s) |29|, stim- 
ulated by the previous similar studies of Ref. |34| , gives less definite estimates than the 
results of applications of our methods (compare the estimate 5R(s) = (— 49^40 ) fl4 HH 



with the result 5R(s) = —97a 4 from Table 2 and Eq. (33)). This fact is related with 
the problems of applicability of the Pade resummation technique to the sign variation 
perturbative series, obtained from the explicit NNLO approximation of the D-function 
with / = 5 numbers of flavours (see also Ref. 0). 



6 Acknowledgements 

We wish to thank K.H. Becks, D. Yu. Bardin and S. Bethke for stimulating the presentation 
of the results of this work prior to publication, at the III Workshop on Artificial Intelligence 
and Expert Systems in HEP (Oberammergau, October 1993), LEP Meetings at CERN (March 
1994) and the University of Aachen (April 1994) respectively. 

We are grateful to R.N. Faustov for attracting our attention to the detailed consideration 



of the results of Ref. |18| at the preliminary stage of our similar QED studies, which will be 



presented elsewhere and to F. Le Diberder for useful comments. 

The work of one of us (V.V.S.) was supported during his stay in Moscow by the Russian 
Fund for Fundamental Research, Grant No. 93-02-14428. 

Note added. Our estimate for d 3 was recently supported by the phenomenological analysis 
of the ALEPH data for R T |||. 



i A similar, more conservative estimate gives 5R T = (78 ± 50)a 4 (for a review, see Ref. |ft3[ ) . 



8 



References 

[1] S.G. Gorishny, A.L. Kataev and S.A. Larin, Phys. Lett. B259 (1991) 144; Pisma ZhetfbS 
(1991) 121. 

[2] L.R. Surguladze and M.A. Samuel, Phys. Rev. Lett. 66 (1991) 560; 66 (1991) 2416 (Erra- 
tum). 

[3] E. Braaten, S. Narison and A. Pich, Nucl. Phys. B373 (1992) 581. 
[4] E. Braaten, Phys. Rev. Lett. 71 (1993) 1316. 

[5] S.A. Larin, F.V. Tkachov and J.A.M. Vermaseren, Phys. Rev. Lett. 66 (1991) 862. 

[6] S.A. Larin and J.A.M. Vermaseren, Phys. Lett. B259 (1991) 345. 

[7] K.G. Chetyrkin, A.L. Kataev and F.V. Tkachov, Phys. Lett. 85B (1979) 277; 
M. Dine and J. Sapirstein, Phys. Rev. Lett. 43 (1979) 668; 
W. Celmaster and R. Gonsalves, Phys. Rev. Lett. 44 (1980) 560. 

[8] B.A. Kniehl and J.H. Kiihn, Nucl. Phys. B329 (1990) 547. 

[9] K.G. Chetyrkin, S.G. Gorishny, S.A. Larin and F.V. Tkachov, Phys. Lett. 137B (1984) 
230. 

[10] S.G. Gorishny and S.A. Larin, Phys. Lett. B172 (1986) 109; 
E.B. Zijlstra and W. van Neerven, Phys. Lett. B297 (1992) 377. 

[11] G.B. West, Phys. Rev. Lett. 67 (1991) 67; Phys. Rev. Lett. 67 (1991) 372 (Erratum). 

[12] V.I. Zakharov, Nucl.Phys. B385 (1992) 452; 

L.S. Brown, L.G. Yaffe and C. Zhai, Phys. Rev. D46 (1992) 4712. 

[13] J. Chyla, J. Fischer and P. Kolaf, Phys. Rev. D47 (1993) 2578. 

[14] J. Chyla, A.L. Kataev and S.A. Larin, Phys. Lett. B267 (1991) 269. 

[15] F. Le Diberder and A. Pich, Phys. Lett. B286 (1992) 147. 

[16] A.A. Pivovarov, Z. Phys. C53 (1992) 461. 

[17] P.M. Stevenson, Phys. Rev. D23 (1981) 2916. 

[18] J. Kubo and S. Sakakibara, Z. Phys. C14 (1982) 345. 

[19] V.V. Starshenko and R.N. Faustov, JINR Rapid Communications 7 (1985) 39. 

[20] G. Grunberg, Phys. Lett. B221 (1980) 70; Phys.Rev. D29 (1984) 2315. 

[21] A. Dhar and V. Gupta, Phys.Rev. D29 (1984) 2822; 

V. Gupta, D.V. Shirkov and O.V. Tarasov, Int. J. Mod. Phys. A6 (1991) 3381. 

9 



J.H. Field, Ann. Phys. 226 (1993) 209. 
R. Barbieri and E. Remiddi, Phys. Lett. 57B (1975) 273. 
M.R.Pennington, Phys. Rev. D26 (1986) 2048. 

S. Brodsky, G.P. Lepage and P.B. Mackenzie, Phys. Rev. D28 (1983) 228; 
Hung Jung Lu, Phys. Rev. D45 (1992) 1217. 

O.V. Tarasov, A. A. Vladimirov and A.Yu. Zharkov, Phys. Lett. 93B (1980) 429; 
S.A. Larin and J.A.M. Vermaseren, Phys. Lett. B303 (1993) 334. 

J.D. Bjorken, report SLAC-PUB-5103 (1989). 

M.A. Samuel, G. Li and E. Steinfelds, Phys. Lett. B323 (1994) 188. 
M.A.Samuel and G.Li, preprint SLAC-PUB-6379 (1993). 
A.L. Kataev, Phys. Lett. B287 (1992) 209. 
K.G. Chetyrkin and A. Kwiatkowski, Z. Phys. C59 (1993) 525. 
P. Nason and M. Porrati, preprint CERN-TH.6787/93 (1993); 

LI. Balitsky, M. Beneke and V.M. Braun, preprint MPI-Ph/93-62, PSU/TH/130 (1993). 
S. Narison, preprint CERN-TH. 71 88/94, PM 94/08 (1994). 

J. Fleischer, M. Pindor, P.A. Raczka and R. Raczka, Proceedings of the XII Warsaw Sym- 
posium on Elementary Particle Physics, Frontiers in Particle Physics, Kazimierz, Poland, 
May- June 1989. Eds. Z. Adjuk, S. Pokorski and A.K. Wroblewski (World Scientific, Sin- 
gapore 1990) p. 412. 

[35] F. Le Diberder, talk at the QCD-94 Workshop, July 1994, Montpellier, France, to appear 
in Nucl. Phys. Proc. Suppl. ed. S. Narison. 



10 



f 

J 


df 


/ Jest\ 

I Uo ) FT'fJ 

\ A J HiKy n 




d% st 


1 


14.11 


7.54 


7.70 


75.4 


2 


10.16 


6.57 


7.55 


49.76 


3 


6.37 


5.61 


6.40 


27.46 


4 


2.76 


4.68 


5.27 


8.37 


5 


-0.69 


3.77 


4.16 


-7.7 


6 


-3.96 


2.88 


3.08 


-20.89 



Table 1 : The results of estimates of the NNLO and NANNLO corrections in the series for the 
D-functions. 



/ 




(Oech 


(Opms 


r .est 


1 


-7.84 


-14.41 


-13.16 


-166.4 


2 


-9.04 


-12.63 


-11.65 


-146.6 


3 


-10.27 


-11.03 


-10.23 


-128.4 


4 


-11.52 


-9.58 


-8.98 


-111.8 


5 


-12.76 


-8.29 


-7.9 


-96.8 


6 


-14.01 


-7.17 


-6.97 


-83.3 



Table 2: The results of estimates of the NNLO and NANNLO corrections in the series for R(s). 



f 


d% £ 






d% st 


1 


44.97 


39.62 


40.78 


423.77 


2 


36.19 


33.28 


34.26 


302.7 


3 


27.89 


27.37 


28.16 


199.61 


4 


20.07 


21.91 


22.50 


113.68 


5 


12.72 


16.91 


17.30 


44.1 


6 


5.85 


12.39 


12.59 


-9.89 



Table 3: The results of estimates of the NNLO and NANNLO corrections in the series for BjnSR. 



/ 


d e 2 x 






d% st 


1 


34.01 


27.25 


28.41 


290.27 


2 


26.93 


23.11 


24.09 


203.37 


3 


20.21 


19.22 


20.01 


129.96 


4 


13.84 


15.57 


16.16 


68.07 


5 


7.83 


12.19 


12.59 


17.79 


6 


2.17 


9.08 


9.29 


-21.57 



Table 4: The results of estimates of the NNLO and NANNLO corrections in the series for 
BjpSR and GLSSR. 
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